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Abstract 

We consider a model for a one-dimensional photonic crystal formed by 
a succession of Kerr-type equidistant spaceless interfaces immersed in a 
linear medium. We calculate the band structure and reflectance of this 
structure as a function of the incident wave intensity, and find two main 
behaviors: the appearance of prohibited bands, and the separation and 
narrowing of these bands. A system with these features is obtained by 
alternating very thin slabs of a soft matter material with thicker solid 
films, which can be used to design a device to control light propagation 
for specific wavelength intervals and light intensities. 
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Photonic crystals PCs — spatially periodic, composite materials that can ex- 
hibit photonic band PB gaps for light propagation — have reached a mature state 
of development, with their optical properties well understood and with many 
realized and potential applications [T]. Most of this research deals with PCs 
whose characteristics are fixed, that is, once they have been fabricated there is 
no possibility to alter their optical response. A recent trend, however, concerns 
tunable or active PCs; by this we imply that, by means of some external agent, 
it becomes feasible to change the optical properties of the PC continuously and 
reversibly. This could lead to tunable optical waveguides, switches, limiters, 
and polarizers; to reconfigurable optical networks; and to electrooptic intercon- 
nects in microelectronics. We can classify tunable PCs according to two broad 
categories. For one of these, an external agent causes structural changes with 
no alteration of the dielectric constants of the constituent materials. In the 
other category the configuration of the PC remains the same, and it is some 
material property of the PC that is affected by the external agent. Structural 
tuning has been proposed or accomplished by means of mechanical stress ap- 
plied to a polymer opal of nanoscale spheres, by applying an electric field to 
a PC on a piezoelectric substrate; by incorporating diodes in a 2D structure 
of wires; and by applying a magnetic field to periodically distributed magnetic 
particles [5] . On the other hand, tuning through the alteration of some material 
property involves the incorporation of a ferromagnetic or ferroelectric material 
in the PC -to be tuned by an external magnetic or electric field, respectively 
[3]. In particular, LCs are well-established electro-optic materials that can be 
tuned by means of pressure, heat, and applied electric or magnetic field. Their 
incorporation within a PC is of particular interest because of the possibility of 
selectively tuning PB gaps, as reviewed recently [4]. 

In this letter we consider a system whose schematic plot is shown in fig. 1. This 
PC consists of a succession of spaceless nonlinear interfaces immersed in a linear 
medium. The nonlinear interfaces are made of a soft matter medium whose 
nonlinear response is that of a Kerr material for which the nonlinear response 
is proportional to the square of the magnitude of the intensity of the light [5]. 
A system with these features can be managed by alternating very thin slabs of 
a soft matter material like a nematic LC, with thicker solid films. Hence, for 
a low power laser beam (for instance He-Ne) the nonlinear index refraction of 
the solid stacks is to be negligible in comparison with that of the nematic films 
whose usual value is more than eight orders of magnitude larger than that of 
a solid |6]. In this model we shall ignore the thickness of the nonlinear slabs 
located in x = na and approximate their optical response by using the Kerr 
model. This allow us to write the following propagation equation 

E Six-na)\A^ = (1) 

\ n— — oo / 

where n2 is the nonlinear refraction index and ip is the amplitude of the prop- 
agating wave. For points x ^ na, the general solution of this equation is of the 
form: V = Ae'^''^ + Be^'^'^^ . However, after a standard integration of Eq.(IT]) 
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around the points x = na [7], we can establish the boundary conditions 

^ (0) = e'^V (a) , e^-^^^ - e^"'^^ = ^^^9^^ (0) IV- (O)l' , (2) 

where we have already introduced the Bloch parameter J5] xa to characterize 
the momentum in each unit cell. Upon applying these conditions to the general 
solution, permit us to write the following transcendental equation 

,. / N cos iz - xa) sin^ z 1 , , , m „ /„n 

f{z,p, xa) = P h COS z - xa) 1 - cos (z + xa ==0, (3) 

z 2 

where z = fca, p = B^gka and 5 = In Fig. 2 we have illustrated the for- 

bidden and permitted bands that we have found by solving Eq.Q numerically, 
where / is plotted in the vertical axis and z is shown in the horizontal one. 

Each plot is the result of the superposition of the different curves obtained 
by varying the parameter xa in the function f (z,p,xa)- Hence, the roots of 
Eq.Q correspond to the points where curves cross the horizontal axis and thus 
each crossing give us a point of the permitted band. Conversely, there is a 
prohibited band whenever there is a white interval (no crossing point). 

Fig. 3 depicts prohibited bands in the interval < 2 < 10 and show how 
they appear and disappear versus the parameter p. 

Some interesting remarks of these plots are the following. Fig. 3a for the 
interval < p < 1 shows that there are no prohibited bands for p = whereas a 
prohibited band appears sd p — 0.03 that widens, until a very narrow prohibited 
band appears in p = 0.3. For the interval 0.3 < p < 0.7, the first band is 
narrowed and the second one is widened. Simultaneously, the first band is 
displaced to the left until the adjacent allowed band located in a; = disappears. 
For 0.7 < p < 0.9, a third forbidden band is added and it also widens against p. 
Fig. 3b for the interval 1 < p < 1000, a fourth prohibited band is added, and all 
of them become narrower. For p — 100 the left prohibited band is split into two 
narrow bands, and for p > 200 there are no more bands. In general terms we 
can address two characteristic behaviors of these bands, one for the appearance 
of the prohibited bands, and another for the separation and narrowing of the 
mentioned bands. 

Instead of an infinite array, for experimental purposes it is wealth to con- 
sider a finite array of nonlinear interfaces n = 1, 2, ... immersed also in a linear 
medium. It is convenient to write the general solution for a nonlinear material 
in the following form 

(x) = c„e'^-^ + d„e-^'=% (4) 

where the subscript n corresponds the stack number n. The boundary conditions 
given by Eq.Q can be expressed in the n-th interface as 

^ ^—ikx I 7 ikx „ ^~ikx i j Akx frr\ 

Cnt + dnC = c„_ie + dn-ie (5j 

and 
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An additional condition stems from the fact that we shall assume that a wave 
travels from the left, and no wave is coming from the right 

do^dn{0) = A, co = Cn{0)^B, c„(Ar)-0 = Cw. (7) 

We solve the system of recurrence equations defined by Eqs.(IS|) and © 
with the initial conditions Eq. ([7]) by straightforward substitution, and calculate 
directly the reflectance defined by i? = \B\ / \A\ . Here we denote by n the 
number of non linear interfaces and iV — 1 the total number of them contained 
in the finite array of stack, in this way the system will have N linear regions. 

Fig. 4a shows the reflectance, R {A, k) as function of the independent vari- 
ables A and k. Fig. 4b exhibits the same information but depicted in a contour 
plot. We present first the simplest case for which n = 1 {N = 2), and g = 0.2 
to get some useful insight. It can be observed that R{A, k) goes to 1, provided 
that A tends to 4, for any value of each k, except for the interval 2.1 < k < 2.3, 
where the material becomes transparent. 

This region disappears when g tends to 1, as one can see in Fig. 5, to reach 
finally the opacity (R {A, k) ~ 1) for g = 6, with = 2. 

Fig. 6 displays a panel plots where g and N are varying, notice that by 
keeping g constant and increasing iV, the trough in the middle of the surface 
is maintained. 

Further calculations allow us to address the following interesting remarks. 
For the case of a only one nonlinear interface {N — 2), R{A,k) increases by 
enlarging A and keeping constant the value of k. Nevertheless, there exists a 
transparency stripe situated approximately in the /c— interval [2,3] which nar- 
rows when g grows until it disappears for 5 = 6. Similarly, for the case of three 
nonlinear interfaces {N — 4): R{A,k) grows for k fixed by enlarging A. Also, 
there is a transparency stripe which is narrowing but it does not disappear 
even for g = Q. The behavior is similar even when the system consists of five 
nonlinear interface (A^ = 6). A general feature observed in all these calculations 
is that by increasing both g and A^, the border of the transparency band is no 
longer smooth. 

We have elaborated a model for a simple nonlinear one-dimensional PC for 
which the nonlinearity is concentrated in certain spaceless barriers periodically 
located. This system can be constructed by piling up slabs of a nonlinear soft 
matter material like a LC inserted between thicker solid films. We have calcu- 
lated analytically the band structure of this periodic structure. We have shown 
that by increasing the intensity of the normal incident wave, there appear more 
prohibited bands which separate each other and get narrow. The same qualita- 
tive behavior is shown by the reflectance of a finite array of alternating stacks 
for which forbidden bands displace and narrow versus the signal intensity. 

Our results suggest to design devices based on this optical structure in order 
to prevent the propagation of light whose intensity is either larger or lower than 
certain threshold, for a specific given interval of the wavelength spectrum. 
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Fig. 1 Schematic plot of the nonhnear photonic crystal. The interfaces exhibit 
Kerr-type optical response whereas the rest of the material displays a 
linear response. 

Fig. 2 Band structure of the nonlinear PC. f {z,p,xo-) versus z parametrized 
by xa. The roots of Eq.Q correspond to the points where curves cut the 
z— axis and thus each crossing give us a point of the permitted band, a) 
p = 0.03, b) p = 0.3, c) p = 0.7 and d) p = 0.9. 

Fig. 3 The clearer intervals denote the forbidden bands versus light intensity p 
(vertical axis) in the interval < z < 10. a) < p < 1 and b) < p < 
1000. 

Fig 4 a) Reflectance R against A and k for n = 1 (A'' = 2), and g = 0.2. b) 
Contour plot of R versus A and k and the same parameters. The clearer 
regions correspond to larger values of R. 

Fig 5 The same as Fig. 4 but for g = I and n = 1{N = 2). 

Fig. 6 Table of contour plots of R against A and k for the shown values of g 
and N. 
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